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Abstract 

We show how perturbation theory may be reorganized to give splitting functions 
which include order by order convergent sums of all leading logarithms of x. This 
gives a leading twist evolution equation for parton distributions which sums all leading 

logarithms of x and Q'^, allowing stable perturbative evolution down to arbitrarily small 
values of x. Perturbative evolution then generates the double scaling rise of F2 observed 

at HERA, while in the formal limit a; ^ at fixed the Lipatov x^'^ behaviour is 
eventually reproduced. We are thus able to explain why leading order perturbation 
theory works so well in the HERA region. 
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The evolution of structure functions at small x is expected to be problematic within 
the conventional framework of perturbative QCD, in that the usual perturbative evo- 
lution appears to be unstable at small x. This expectation rests on the observation that 
singlet anomalous dimensions grow in an unbounded way as — (specifically, at 
order m in perturbation theory 7^(a) ~ (;^)"^)- Indeed, the leading order term drives 
a corresponding growth in the singlet structure function as both ^ and increase 0. 
However if ^ increases much more rapidly than t = InQ^/A^, higher order terms be- 
come increasingly important, and it is no longer clear whether perturbation theory can be 
trusted [H- Hence it is often claimed that the usual evolution equations |jl|] must then be 
abandoned in favour of equations which resum the logs of ^ at fixed t, resulting in an 
even stronger (power-like) growth. 

Strong scaling violations at large accompanied by a corresponding increase at 
small values of x are indeed observed in recent measurements of the proton structure at 
HERAp,^. The leading order pertiirbative prediction 0] can thus now be tested empir- 
ically: the structure function data should display a double scaling behavior [0 in the two 



variables a ~ yln^lnt and p ~ yln^/lnt, rising linearly in the former while remain- 
ing relatively independent of the latter. This scaling follows directly from the fact that 
the leading order small-x evolution equations reduce asymptotically to an isotropic wave 
equation in the plane of In ^ and Int. In fact not only is the scaling observed 0, but even 
the slope of the rise turns out to be remarkably close to the leading order prediction 0. 
These results remain stable upon the inclusion of two loop corrections , and it now seems 
possible to detect two loop effects in the data through a small reduction in the slope. 

All this suggests that it should be possible to keep the behaviour of the perturbative 
expansion under control by summing leading logs of - as well as logs of , while remaining 



within the standard framework of perturbative evolution [|l0[ . The double asymptotic scal- 
ing behaviour observed in the HERA data should then emerge as the dominant behaviour 
in some well determined region of the x-t plane. Such a programme is possible since the 
leading singularities in both and 7^ have now been explicitly computed to all orders 
in Us/N [^T[|T^. However in a previous study |]T^, in which three and four loop singulari- 
ties were retained, it appeared that instabilities may develop which signal the breakdown 
of perturbation theory. Here we will explain in detail how the perturbative expansion of 
anomalous dimensions can be reorganized at small-x in a way which is consistent with 
the renormalization group, but in which the leading logs of ^ are summed to all orders. 
We will then prove explicitly that this leads to stable perturbative evolution equations for 
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parton distributions at all positive values of a;, in which all corrections are either of higher 
order in as or are higher twist, and which thus breaks down only at low values of Q^. 
Solution of these equations will then enable us to determine in precisely which region of 
the x-t plane double scaling should remain valid, and whether it will be possible to observe 
the stronger power-like rise characteristic of solutions of the BFKL equation. 

According to the operator product expansion and renormalization group, the Mellin 
moments /^(t) = Jq dxx^ f'^{x; t) of the parton distribution functions t) evolve mul- 
tiplicatively with t: 

j/Nit) = Y.^%{as{t))f^{t), (1) 
j 

where 7]v(ct) are the anomalous dimensions of local operators 0]. The perturbative evolu- 
tion then sums all logarithms of Q^, so the 0{a) and 0{a^) terms are known as the leading 
and subleading log approximations respectively. Perturbative expansion of the anomalous 
dimensions is justified by factorization theorems; these have been shown recently to apply 
to all orders even at small x [|T3| - p!5l . The inverse Mellin transform of gives 

|r (x; ^) = E ^yp^'^y^ t)f\-y-^t). (2) 

where the splitting functions P^^{x;t) are related to the anomalous dimensions by0 

l%Mt))= fdxx''P'\x-t). (3) 



This version of the evolution equation is more physical because it only involves integration 
over the physically accessible region y > x . 

At small X it is of course necessary to take into account the fact that, if ^ is large 
enough, Feynman diagrams which contain powers of In ^ may be just as important as 
those with powers of InQ^. It is thus no longer appropriate to organize the perturbative 
expansion of the splitting functions in powers of as only. In the Mellin space equation 
(P these extra logarithms are due to the presence of poles in the anomalous dimensions 
as ^ 0. Now, while nonsinglet anomalous dimensions are regular as A ^ 0, singlet 
anomalous dimensions computed at leading order in a have a singularity; aX m—th order 
they have at most (at least in reasonable renormalization schemes) an A""* singularity 



^ This definition difi^ers from the usual one by a factor of also our moment variable A'^ 
differs by one from the usual one so that for us the first moment is that with N = 0. 
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pl] - |T5|1 . Hence, separating out these singularities, and suppressing the flavor indices i,j 
for clarity, the anomalous dimensions admit the perturbative expansion (see fig. |l|a) 



oo m oo ^ m \ 

m=l n= — oo m=l ^ 71=1 ' 



(4) 



where are simply numerical coefficients, and 7^"* are all regular as — > 0. This then 
corresponds to the expansion 

Pir, t) = ± {car E + ^'""(^)) 

m=l ^ n=l ^ ' 



of the splitting function, where P^'^\x) are regular as a: 0. At small x the usefulness 
of this expansion in powers of a is spoilt by the logs of 1/x which can compensate for the 
smallness of as{t). 

However this is not the only way to order the expansion: if instead of expanding in 
a and then in A^, as in (^), we choose instead to expand in a and a/N (see fig. |lp3), the 
anomalous dimensions take the form 

oo oo / \ ^ 

m=l n=2—m ' 

which corresponds to the splitting function 

[n-V]\ ~ 



p(.;.)^f:M0)"'-'(iE^r.--'^^^if3^?; 

m=l n=l 

+ i;^-r-'Kw)"'^«(i--)) 



(7) 



While not so useful at large x, at small x this is clearly the most appropriate expansion, 
since at each order in (Xg all the leading logs of 1/a; have been summed up. 

Solving the evolution equations with the usual expansion eq. (|5|) of the splitting func- 
tions corresponds at leading order (m = 1) to summing up all logs of the form 

alihiQ'yilJ^ (8) 

with q = and < r < p. If instead only the sum of leading singularities eq. (^ (again 
with m = 1) is included in the splitting functions, then solving the evolution equations 
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sums the leading logs (|^) with r = p, and 1 < q < p- Of course, if contributions which 
are higher order in as are included in either expansion of the splitting function, eventually 
all leading logs of both ^ and are included. However while in the former logs of 
are considered leading, in the latter logs of ^ are leading, the roles of q and r in (||) 
being interchanged. Both expansions are consistent with the renormalization group in the 
sense that a change in the scale at which the m-th order contribution to the expansion is 
evaluated is equivalent to a change in the m + 1-th order term. 

It is now easy to see how to construct an intermediate "double leading" expansion 
which includes all leading logs, i.e. one such that perturbative evolution sums all terms 
with l<q<p,0<r<p, and 1 < p < q+r, so that each extra power of as is accompanied 
by a log of either ^ or or both. To do this, the anomalous dimensions are expanded as 
(see fig. |l]c) 

CO / ca / \ \ 

m=l ^ n=l ^ ^ ' 

Each subsequent order of the expansion contains an extra power of a in comparison to the 
previous one, so the scheme is still consistent with the renormalization group. Of course, 
any number of renormalization group consistent expansions which interpolate between the 
double leading one (eq. (^) and the large x (eq. (^) and small x (eq. (^)) ones can also be 
constructed. The important point is that if we choose an expansion which is appropriate at 
small X (say, the extreme small x one (§[)), then each subsequent order of the expansion is 
genuinely of order ol as compared to the previous one, all logarithms having been included 
in the coefficients. When the scale increases, the higher order contributions are then 
asymptotically small. 

The double scaling behaviour described in refs. [^^, and observed at HERAp,^, is 
governed by the pivotal m = n = 1 term in the various expansions (§),(^ and (|^), with 
(A]^)^^ = 7^/3o/47r determining the slope of the rise of at small x and large . The 
subleading terms determine b (as defined in ref. 0) and the relative normalization 
of the quark and gluon distributions, while A| = 0: the leading two-loop term A\ is 
responsible for a small but observable reduction in the steepness of the slope Thus at 
the intermediate range of x and t currently being explored at HERA the double leading 
expansion (^) is the most appropriate, and should be used for a detailed comparison to 
the data. Here we will instead be interested in the small x limit and will thus henceforth 
concentrate on the small x expansions eqns. (^) and (^, which are sufficient to demonstrate 
the novel features of perturbative evolution in this kinematic regime. 
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Before proceeding further, we must confront the convergence issue with which we 
began the paper. Consider first the leading term in the series (^): since both 7^ and 7^ 
vanish at leading order, audi 7^ = (Cf/C^)7^ we consider explicitly 

00 

7^^(«) = E(^n)^^ {^y+0{a) ^ {4ln2)-'A{X/N)+0{a), (10) 

n=l 

where 

A(a)=41n2f^Y (11) 



TT 



Now, it is well known (see e.g. ref. [^) that this series develops a branch-point singularity 
at = A [3, so that while the series converges for |A^| > A, it is not even Borel summable if 
Re N < X . This seems to pose an insurmountable problem for the perturbative approach 
to small X evolution: the series which defines the leading coefficient in expansion (|^) of the 
anomalous dimension, which was supposed to be useful for small A^, is actually only well 
defined when is sufficiently large. This apparent inconsistency seems to have led many 
to the conclusion that conventional perturbative evolution breaks down at small x. 

In fact the dilemma is resolved by the observation that the physically relevant quantity, 
namely the splitting function P^^(x), is instead given at leading order in the expansion 
(^ by the series 

X ^-^ (n — 1): -K X 

n=l 

where now As(t) = A(as(t)). This series converges uniformly on any finite intervals of x 
and t which exclude x = 0. To prove this, we write the two series (^0|) and (|T2|) as 

cx) 00 

A{v) = J2 <^nv\ B{u) = J2 r^^""'' (13) 

where the coefficients a„ = 41n2(A^)^^(41n2 Ca/tt)""^. The series for A(v) then has 
radius of convergence one, from which it follows trivially that the radius of convergence of 
B{u) is infinite. 

The defining relation (0) may now be rewritten in the form 



POO 

A{v) = / due-''/'' B{u) 
Jo 



(14) 



^ Following the usual conventions, Ca = 3, Cf = |, Tr = ^ for QCD with three colors. 
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B is the Borel transform of A. The series A{v) is not Borel summable for Kev~^ < 1 (i.e. 
ReA'" < A) because the integral no longer converges at the upper limit (corresponding to 
the lower limit of the integral over x in (^)): the real coefficients are all positive and 
decrease only very slowly as n — > oo, so B{u) ~ as 1^1 oo.i It follows that the only 
reason for the bad behaviour of the series ([1^) is that when transforming to Mellin space 
one attempts to integrate all the way down to a; = 0, and this is not possible for singlet 
distributions because the total number of partons diverges there. If instead the parton 
distributions are evolved using the Altarelli-Parisi equations (H), the splitting functions 
are only required over the physically accessible region x > x^^i^ and no convergence 
problems arise. Indeed for all physical applications it is sufficient to truncate the series 
( IT2I) after a finite number of terms .@ 

Similar considerations presumably apply to subleading terms in the series (^ and (|^). 
Recently the first nonvanishing term in the large expansions of 7'^^ and 7'^'^ have been 



calculated [|I4[: the result may be written in the form 



7^^(«) = Tr^{1 + ^(A/iV)) + 0(a), 7^^(«) = ^^^M>^/N) + 0{a), (15) 

where A{y) has the same form as A{y) (p!3|), with new coefficients dn, but the same radius 
of convergence. It follows immediately that the corresponding expressions for P'^^{x;t) 
and P'^'^{x;t) are convergent for all nonzero x. It is tempting to conjecture that the same 
must be true for all the coefficients in the small x expansion (0): aU that IS necessary 
is that the coefficients in the expansion of the corresponding anomalous dimensions 
have a positive radius of convergence, i.e. that there is no singularity (or accumulation of 
singularities) at = 00, or equivalently at a = 0. This is a conventional assumption in 
perturbative QCD. 



^ For timelike anomalous dimensions, relevant for the small x evolution of fragmentation func- 
tions, the corresponding series have alternating coefficients, and thus converge throughout the 
physical region. The same is true of (pg theory. 

^ In practice this means that when solving the evolution equations (|l]) numerically, the di- 
vergence of the series expansion of the anomalous dimensions at small N can simply be ignored: 
although adding another term in the series ([To| ) can change the anomalous dimension a lot at small 
A'', this will be compensated by a corresponding shift in the steepest descent contour, such that 
the change in the evolved distribution is actually very small. We have checked this numerically. 
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Using the analytic expressions given in [|TT,n] we have computed the coefficients an 
and dn to very high orders. In practice only the first dozen or so are needed in any realistic 
calculation: we give the first thirty six in the table. 

We are now finally able to write down evolution equations for parton distribution in 
the leading order of the small x expansion of the splitting functions . Because only the 
P^^ and P^'^ splitting functions are nonvanishing at leading order, only the Altarelli-Parisi 
equation satisfied by the gluon distribution survives; the determination of the quark dis- 
tribution, and thus of structure functions, requires the inclusion of subleading corrections 
and will be discussed at the end of the paper. As in the evolution equation may be 
simplified by introducing the variables 



and defining G'(^, () = xg{x; t). Using the splitting function (|T2|) in the evolution equation 
(ID then gives 



9C Po (n-l)! J_^^ 

2 



Y^an+iXsiCr d^2... dC'GiC'X)- 

, n "^—^0 J —^0 J—^0 



[17) 



7 

n=0 



Here 7 = ^/ACaJPo (as in [0)- have used the one loop form of as{t) = Po 



11 — 3^-/, as is appropriate for a leading order calculation. Since we retain only singular 
contributions to the splitting functions the lower limit ^0 = In — in the integrations on 

Xq 

the right hand side can be consistently set to zero. 

Note that if the large ^ (i.e. small x) limit is approached by letting ^ grow as ('^ with 
any /c > 1, but more slowly than e'', corrections to the leading order splitting function are 
exponentially suppressed. If, however, we let ^ grow faster than e'' (i.e. 1/x grow faster 
than a power of Q^), then higher orders in the leading singularity expansion of the splitting 
function will eventually become significant, growing as e^^" 

Differentiating both sides with respect to ^, the gluon evolution equation eq. (|l^ can 
be cast in the form of a wave equation as in ref. : 

fl2 °o /-Si 



d^dC 



G{C, C) = 7'G(e, + 7' E ^n+iXs{Cr d^i d^2... d^' G{C', 

Jo Jo Jo 



(18) 



-I'Gi^.o+i'f^an+i^^^ I rfe'(e-er-^G(e',c). 
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If only the first leading singularity is retained, then eq. (^8|) is simply the wave equation 
discussed in ref. (with 6 = 0; subleading corrections have not yet been included), and 
leads directly to a double scaling behaviour of the gluon distribution when the boundary 
conditions are sufficiently soft. When the higher singularities are also included on the 
right hand side, then the generic features of the propagation in the plane due to the 
wave-like nature of the evolution are preserved. In particular, the propagation is still linear 
and causal, and far from the boundaries independent of the detailed form the boundary 
condition. 

For any reasonable values of ^ it is clear that only the first few terms in the expansion 
of eq. ( IT2D will contribute (notice that in fact 02 = as = as = 0) , and then the asymptotic 
form of the solution will be essentially unchanged (i.e., it will display double asymptotic 
scaling). The easiest way of seeing this explicitly is to recall that the leading order solution 
in which only the first term on the right hand side is retained is just given by a Bessel 
function G'o(^, C) ~ -^0(270"), where, as in 0, we find it useful to define the scaling variables 
cT = p = a/^/C- The asymptotic behaviour Io{z) ~ z~^l'^e^ as z — 00 then gives rise 
to double asymptotic scaling. If this is substituted back into the new terms on the right 
hand side, to generate the first term in an iterative solution, we may use the fact that 
z'^ln-\{z) = -^{z'^In{z^ to writc ([T8| ) in the approximate form 



92 



GiC, ^ 7'G(^, + 7' E ^-+1 (^^) ^-(27^). (19) 



Since p^In{2'ya) is bounded above by ^"^/o(27(j) it follows that double scaling will always 
set in asymptotically provided the series "^n+i (^^fjf^)"^ converges uniformly, that is 

provided ^ < x^c) ■ 

However for very large ^ the splitting function will eventually be dominated by the 
higher orders of the series eq. ([T^). This is the limit which is often approached by means 
of the BFKL equation which resums all leading logs of ^, i.e. all terms with p = r, 
g = in (^), but with the dependence of suppressed. Since at leading order eq.([T^) 
already includes all the leading twist information contained in the BFKL equation (but 
with all higher twist and infrared behaviour factored out), it should reproduce the same 
power-like rise in the limit ^ ^ 00, provided that we freeze the coupling. We can both 
check this, and understand how the behaviour of the gluon evolution changes when the 
coupling runs, by taking advantage of the fact that when the series eq. (|T2|) is dominated 
by its higher order terms, the evolution equation eq. (ffSl) takes a simple closed form. 



Because the series (|T0|),(^3|) has unit radius of convergence it foUows that hmn oo 
!t^±l. — > l.i But then setting On+i ~ in the sum in eq. (ffSl), shifting the summation 
index by one unit, and using eq. (|T7|), we have 



n=l 



^0 ^0 



i n — 1 



dG 



72^a.+iA.(C)^ I d^i I d^2... I d^'G{^',0 = UO^- (20) 



Hence, asymptoticaUy as ^ ^ oo the evolution equation (|T^) becomes simply 

^ dC 



Q^Q^Gi^X)-Xs{0^ = 7'G{CX). (21) 

and asymptotically the summation of all leading singularities leads to a damping term in 
the wave equation.i 

It is now trivial to recover the singular 'Lipatov pomeron' behaviour: if the coupling 
is frozen, then As is just a constant. A, and the solution to eq. (^T]) is given in terms of the 
(double scaling) solution G'o(^,C) to the original wave equation 

Git = e^^Goit = x-^Goit 0- (22) 
The solution with fixed coupling thus displays a strong power-like growth Q in the region 

2 

^ ^ ^C, while below this line we have the usual scaling behaviour (fig. |]a). This situa- 
tion is very similar to that obtained by solving the original A = equation with a hard 
boundary condition ^(^,0) ~ e^^ [p!8|,|7|JTm ; however now the hard behaviour is generated 
perturbatively, by the singularity at = A in the anomalous dimension eq. (|T^), rather 
than imposed by hand. 

Even though (^T]) may be solved exactly for fixed coupling, it will be useful for the 
sequel to derive the asymptotic behaviour of the solution (^2]), using the usual technique 
of Laplace transformation and steepest descent. Defining Q{s, Q = J^d^ e~^^G{t (|21] ) 
may be written as 

^^(«,C) = ^^(«,C) (23) 

which with A fixed integrates to 

\ng{sX)=lng{s,0) + -^C- (24) 

s — A 



^ It can be seen from the table that this asymptotic form begins to set in only for n ~ 20 
however. 
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As was anticipated, but not proven, in ref. Q. 



9 



Estimating the Bromwich integral by steepest descent, the saddle point is located at sq = 
A + 7/p, and thus for a soft boundary condition 0) ~ 1 we find that as ^ ^ 00 

For a hard boundary condition 0) ~ the first factor would be absent: when 
leading singularities are included, the form of the boundary condition is only important 
very close to the boundary (unless of course the boundary condition were even harder than 
x~'^). This feature persists in solutions to the full equation (|T^). 

Of course there is really no justification at all for freezing the coupling: we must solve 
eq. (^1]) with running coupling A(C) = Aqc"^. Since an exact solution is no longer available, 
we work with the with the Laplace transformed equation (^3]) , which now integrates to 

(26) 



In g{s, = In ^(., 0) + ^ (C + In (^^^^) 



and we find a branch cut from X{Q to Aq. The saddle point condition leads to a tran- 
scendental equation, so we now treat the large and small ^ limits separately. When s 
is large, we can ignore the cut, expand out the logarithm, and find a saddle point at 
So = ^ — + • ■ ■, which gives the asymptotic behaviour 

G(e, C) ~ J-e^^'^+^^o-xiOP^ (27) 

^/c^ 



i.e. double scaling up to a small correction. This holds throughout the region where 

2 

the correction is small, i.e. for ^ <^ jx^zryy^C^- The behaviour for large ^ is found by 

dominating the Bromwich integral by the branch point singularity at s = Aq: the saddle 

2 

point is then at sq = Aq + + • • which gives 

^(^,0 - ^(^(Ao - X{0)f^''e^''+^'</'^. (28) 

Again this behaviour is valid whenever the corrections to it are small, which means that 

^ > Thus the running of the coupling leaves the exponent of the power-like growth 

^0 

unchanged (at Aq), but severely limits the region in which this is the dominant behaviour 
(see fig. ^). The region in which the scaling behaviour (^7|) obtains is correspondingly 
increased. 

In practice the damped wave equation (^TJ) is still a rather poor approximation to the 
full equation (^Sj), since, as may be seen from the table, the ratio ^ ^ 1 for all n > 1. 
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A better approximation may thus be obtained by setting = eai; this gives the rather 
more comphcated equation 

which however on Laplace transformation becomes simply 

The asymptotic behaviour of the solution may be found just as above: for running coupling 
( PUD again integrates to (^B|), but with the logarithm now suppressed by a factor of epsilon, 
which means that the double scaling solution (|2^) holds throughout the even larger region 

2 

e2(Ao-A)^^^' subleading term in the exponent being suppressed by e. The power- 
like growth ( PBD is now confined to the region ^ ^ (Xq-x) ^'^^^ ^^^^ H^-'' ^^^^^ ^ ^ O--*-' 
this means that unless ( is very small, the Lipatov growth only sets in when ^ is exceedingly 
large.0 

These results are confirmed by a numerical analysis of the full evolution equation 
Eq. ([1^). The equation was solved by using the integrated form 

G(?,C) = Go(?,C)+ /V /V/o(27V(?-a(C-C')) 

^1 (^-l)Wo 

where G'o(^,C) is the solution of the leading order equation, as given in [0. In practice 
the series on the right hand side converges to an accuracy of less than one per cent after 
only a dozen or so terms. Eqn.(|3T|) may be solved by the standard iterative procedure: 
convergence is achieved in practice after only a few iterations. For definiteness we chose 
72/ = 4, A = 120MeV (so that as{M^) = 0.116), soft boundary conditions G{^,0) = 
G{0, C) = const, and xq = 0.1, Qq = IGeV. The result is shown in the contour plot fig. ||a, 
to be compared with the leading order plot fig. 3a of ref. [|^. Even though the range of ^ 
has been increased, there are no significant scaling violations except very close to the left 
hand axis. 

^ For example if C ~ 0.4 (corresponding to around lOGeV^), we would need ^ to be of order 
100. 
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In summary, our analysis of the small x all-order gluon evolution equation eq. ([Tq ) 
shows that the leading order double scaling behaviour holds asymptotically throughout 
most of the ^-^ plane, the power-like singular behavior eq. being confined to a very 
small wedge close to the boundary C = 0- This is to be contrasted with the situation for 
fixed coupling 0] , or the phenomenological approach in which a singular behaviour is input 
to the leading order evolution equation [lT8|,[7|,p!6| : in both cases the power like singularity 

, 2 

is preserved by the evolution and propagates into a large region ^ ^ 'x^C- The reason that 
the region in which the power-like growth develops is in reality so small and narrow is 
basically twofold: the running of the coupling means that the growth near the boundary 
is rapidly damped away by the exponential fall of ctg as ( increases, and the size of the 
growth is in any case severely limited by the smallness of the coefficients of the leading 
singularities. It thus seems extremely unlikely that the 'hard pomeron' will be seen in 
structure functions, since at any reasonable value of Q^, the power-like rise only sets in 
at extremely small values of x. In ^ { "a^{t) Y^^^ 1/e ^ 20, and in any case such small 



values of x are well inside the recombination region [ff9|,|T5| In - ^ In — ( — ttv)^, Xr ~ 10 
where perturbation theory has already broken down, and the rise is presumably damped 
by the nonperturbative effects necessary to restore unitarity [Q. 

Having found the gluon distribution, we may now discuss the determination of the 
quark distribution and the structure function F2{x]t), given in the parton scheme by 

F|(a:;t) = 4g(x;t) + F2^^(x;t) 

^ (32) 
Q{x;t) =x2_^ {qi{x]t) +qi{x;t)). 

In the small x region, the nonsinglet contribution F^^{x;t) may be neglected. Since 
7^^ and 7'^'^ begin at order a (see (^5])) to calculate Q from G we must go to next-to- 
leading order in a in the expansion eq. (^. Writing this as 7*-' (q;) = J2T=o'^'^'ya ^ 
diagonalize 7*-' (a) order by order in a. To first order the eigenvalues are 
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gq 

A+(a) = 7o^^ + a(7r + ^7^ + 0{a') 

7o 

gq 

7o 



(33) 



The corresponding eigenvectors are given by 



qg gg 

Q+ = a\G+ + 0{a^) Q" = -^G" + 0(a). (34) 
7o 7o 
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The eigenvectors are not orthogonal, because 7*-' is not symmetric. At leading order the 
larger eigenvalue A"*" leads to the wave equation (ffSI) for , while the small eigenvalue 



vanishes, so G and Q do not evolve. At order a, ( |T5D gives A (a) = — ^ '^cj^ 



so 



G-{t) = G]^{0)e-^'<, = Q^(0)e-'^''^, where d' = ^f3oK The Ml correction to A+ 

is as yet unknown, since only the first two terms of have been computed. Retaining 
only the first (one loop) term the evolution equation (^Sj) is supplemented by a damping 
term proportional to 5 = (11 + ^^)/(3o as in 0. The solutions G~^{C) and Q~^{C) then 
acquire an extra factor of e"'^''. 

Using the eigenvalues (^3[) and eigenvectors (0), and the expansion (^5]) for 7^^, it 
is now not difficult to derive the quark evolution equation: taking as boundary conditions 

G{^,0)^Go{O,Q{C,0)^Qo{O, 

^Q(e, = -S'QoiOe-''^ + [g(^, C) 

°° tMn-l _ (35) 

+ J]a.A.(C)- / de\^^G{eX) , 
n=i Jo [n-L). J 

where the gluon distribution G{$,X) is a solution of the wave equation (p!8| ) with boundary 
condition G'(^, 0) = G'o(^) + |Qo(0- This equation may then be readily integrated to 
give the quark distribution and thus -F|'. Although we have derived this equation in the 
parton scheme, it is actually scheme independent: changes in the renormalization scale 
only affect terms which are of higher order in Q!s.i In this sense the expression (^5]) is still 
leading order, and should thus be treated on the same footing as ( [T8|) ; scheme dependent 
subleading corrections first come in through higher loop singularities in 7^^. 

We may now discuss the asymptotic behaviour of the quark distribution. Since 
throughout most of the plane G(^, C) ~ A/'/o(27a), we may use the same trick as 
we did to derive ( p^Qj ) to write ( |35| ) as (setting d and d' to zero for simplicity) 



^Q(?,C) - ^7^Arf;a.(^') /.(27a). (36) 



^ This may be proven explicitly; although in a different scheme (such as MS) the coefficients 
a„ will be modified, the difference is made up by the new 0{as) corrections to the coefficient 



functions [14|. 



13 



It follows that when ^ « ^Ce^^, C) scales just as observed experimentally ffl. For 

2 

^ ^ ttCc^^, if we set = 1 for all n, (^61) may be simplified yet further by using the 
relation J2°^oo ^"-^n(^) = exp {^z{t + to give 

Q(^,C)~A^Ce«"°+'^'^/^°. (37) 

This result may be readily confirmed using Laplace transforms with respect to ^. 

It follows that even when the gluon scales, F2 can still grow rapidly sufficiently close to 
the C = boundary because of higher order singularities in the anomalous dimension 7^. 
Indeed these singularities turn out to be much more important than those of 7^, despite 
being suppressed by an additional power of a, essentially because the coefficients dn ^ an- 
Note that in the BFKL equation quark loops are suppressed, so the contribution of 
these singularities is not included. In fact even if it were possible to detect the tail of a 
power-like rise in the gluon distribution, in the structure function this would be masked 
by a similar rise generated by the subleading singularities in the coupling of the gluons to 
quarks. This effect may be confirmed numerically by evaluating (^5]), using as input the 
gluon distribution displayed in fig. and taking as boundary condition Qq ~ ^Gq. The 
resulting contour plot of In F2 is displayed in fig. ^5. Indeed F2 differs very little from the 
double scaling result except in the wedge close to the left hand boundary, where instead 
the behaviour (Wfh is found. 



The growth (^) should not be taken to literally, however, since for any reasonable 
value of X it only occurs for values of very close to the starting scale Qq, at which 
nonperturbative effects are also relatively important. Furthermore, in this region higher 
twist recombination effects are expected to become significant . More interestingly, the 



higher order terms in increase the relative normalization of Q with respect to G at 
small X for a given starting scale Qo (by around a factor of two). Since when Qq = 1 GeV 
the leading order result has approximately the correct normalization, this suggests that at 
all orders the starting scale Qo at which the soft pomeron boundary condition is imposed 
should be raised to compensate. In fact, a starting scale of around 2 GeV turns out to be 
required; the resulting contour plot of lnF2 is shown in fig. Since Aq is a little lower, 
the scaling behaviour is also improved, the power-like growth being now confined to the 
top left-hand corner of the plot. 

To make the comparison with the double scaling results of ref. more quantitative, 
we may compute the slope of the linear rise of (in the notation of [0,^) lni?^F2 in the 
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HERA region. In leading order this is simply given by 27 = 2.4, which is reduced by around 
20% by two loop corrections 0, whereas the HERA data give 2.4 ±0.2. The corresponding 
slopes for the three distributions plotted in fig. |^ are 2.5, 2.9 and 2.6 respectively. Clearly 
a detailed fit to the data would thus require the inclusion of all of these competing effects. 

In conclusion, the leading twist evolution equation which sums all leading logarithms 
at small x, but is also fully consistent with the renormalization group, interpolating 
smoothly between small x and large x, is just the Altarelli-Parisi equation (^), with 
the splitting function expanded in powers of as{t), but retaining terms to all orders in 
as(t)ln^. Perturbation theory at small x is then placed on the same footing as pertur- 
bation theory at large x: there are no instabilities or problems of convergence. In fact 
deep inelastic scattering at small x and large is a much cleaner perturbative environ- 
ment than more traditional deep inelastic scattering at large x, since the non-perturbative 
component of the structure function is relatively insignificant 0,^. Taken altogether, the 
analysis presented here explains why the simple double asymptotic scaling picture pre- 
sented in 10], though only at leading (and subleading 0) order in a^, accounts so well 
for the rise in F|' seen at HERA [|],^ . Perturbation theory breaks down only when higher 
twist effects become necessary to restore unitarity. 

Although they have little effect on the qualitative form of F2 for any reasonable 
values of Q^, the higher order logs in P'^^ do make a significant impact on the relative 
normalization of the quark and gluon distributions. This means that attempts to extract 
the gluon distribution from F2 using the standard two-loop evolution equation should 
be treated with caution: an all-order calculation could give a gluon distribution which 
is smaller by as much as a factor of two. The summation of subleading logarithms will 
become increasingly important as the quality and range of structure function data at small 
X improves further. 

Acknowledgement: We are very grateful to S. Catani for several useful discussions 
on this subject. 
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Figure Captions 

Fig. 1. The terms summed in tlie various expansions of the anomalous dimensions and 
associated splitting functions: a) the standard (large-a;) expansion (^), b) the 
small- a; expansion (^), and c) the 'double-leading' expansion (^). Leading, sub- 
leading and sub-sub-leading terms are indicated by the solid, dashed and dotted 
lines respectively; m denotes the order in a, while n denotes the order in 
Singular terms are marked as crosses, while terms whose coefficients known at 
present (for 7^) are marked by circles: the term which leads to double scaling is 
marked with a star. 

Fig. 2. Regions of validity in the ^-^ plane of various asymptotic behaviours: a) from 
the damped wave equation ( |2T|) with fixed coupling, b) the same with running 
coupling, and c) from the equation ( [2^ ) which takes into account the smallness of 
the coefficients a^. The region marked 'S' denotes the double scaling region, while 
'L' denotes the Lipatov (power rise) region, and 'R' the recombination region. 

Fig. 3. Contour plots of the logarithms of the solutions of the small- a; evolution equa- 
tions in the plane: a) the gluon distribution G, b) the structure function F2 
with Qo — 1 GeV, and c) the same but with Qo = 2 GeV. The values of the 
other parameters are explained in the text. The contours are all equally spaced: 
the interval between adjacent contours is AlnG = 0.85, AlnF2 = 0.62, 0.58 
respectively. A scatter plot of the HERA data0,|[ is superimposed for reference. 



18 



n 






1 


1 


0.78145981 


2 





0.29913310 


3 





0.38806675 


4 


0.11279729 


0.25256339 


5 





0.17838311 


6 


0.01265760 


0.22632345 


7 


0.03816969 


0.15473331 


8 


0.00160119 


0.13891462 


9 


0.01142194 


0.15744938 


10 


0.01742873 


0.11602071 


11 


0.00260717 


0.11532868 


12 


0.00888439 


0.11997040 


13 


0.00942680 


0.09568020 


14 


0.00300214 


0.09846280 


15 


0.00670932 


0.09707906 


16 


0.00581292 


0.08291446 


17 


0.00301620 


0.08559381 


18 


0.00507784 


0.08203738 


19 


0.00400435 


0.07384063 


20 


0.00283184 


0.07550427 


21 


0.00390386 


0.07157467 


22 


0.00301599 


0.06682914 


23 


0.00256389 


0.06747701 


24 


0.00306948 


0.06392641 


25 


0.00242542 


0.06112241 


26 


0.00227736 


0.06101759 


27 


0.00247526 


0.05808061 


28 


0.00203856 


0.05633377 


29 


0.00200551 


0.05576012 


30 


0.00204715 


0.05343445 


31 


0.00176238 


0.05224401 


32 


0.00176276 


0.05142727 


33 


0.00173286 


0.04961887 


34 


0.00155085 


0.04871349 


35 


0.00155308 


0.04780818 


36 


0.00149656 


0.04640235 



Table: The coefficients an and a„, computed us- 



ing formulae in ref. [11| and ref. [14| respectively. 
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